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We have investigated particle acceleration and shock structure associated with an unmagnetized relativistic jet 
propagating into an unmagnetized plasma. Strong magnetic fields generated in the trailing jet shock lead to 
transverse deflection and acceleration of the electrons. We have self-consistently calculated the radiation from 
the electrons accelerated in the turbulent magnetic fields. We find that the synthetic spectra depend on the 
bulk Lorentz factor of the jet, the jet temperature, and the strength of the magnetic fields generated in the 
shock. We have also begun study of electron acceleration in the strong magnetic fields generated by kinetic shear 
(Kelvin-Helmholtz) instabilities. Our calculated spectra should lead to a better understanding of the complex 
time evolution and/or spectral structure from gamma-ray bursts, relativistic jets, and supernova remnants. 



1. INTRODUCTION 

Recent kinetic simulations have focused on mag- 
netic field generation via electromagnetic plasma 
instabilities in unmagnetized flows without veloc- 
ity shears. Three-dimensional (3D) particle-in-cell 
(PIC) simulations of Weibel turbulence [Nishikawa et 
al. 2005, 2009a] have demonstrated subequipartition 
magnetic field generation. We have calculated, self- 
consistently, the radiation from electrons accelerated 
in the turbulent magnetic fields. We found that the 
synthetic spectra depend on the Lorentz factor of the 
jet, the jet's thermal temperat ure, and the strength of 
the g enerated magnetic fields [Nishika wa et alJl201lL 
l2012j . 

Velocity shears also must be considered when study- 
ing particle acceleration scenarios, since these trigger 
the kinetic Kelvin-Helmholtz instability (KKHI). In 
particular the KKHI has been shown to lead to par- 



ticle acceleration and magnetic field amplification in 
relativistic shear flows [Alves et al. 2012; Liang et 
al. 2012]. Furthermore, a shear flow upstream of 
a shock can lead to density inhomogeneities via the 
MHD Kelvin-Helmholtz instability (KHI) which may 
provide important scattering sites for particle acceler- 
ation. 



2. THE STANDARD RADIATION MODEL 

A synchrotron shock model has been widely 
adopted as describing the radiation mechanism 
thought responsible for observed broad-band GRB af- 
terglows [e.g., Zhang & Meszaros 2004; Piran 2005a; 
Zhang 2007, Nakar 2007]. Due to the lack of a first 
principles theory of collisionless shocks, a purely phe- 
nomenological approach to the model of afterglow ra- 
diation has been prescribed. Firstly, electrons are as- 
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Figure 1: Figure [Tk shows the spectra for 
the cases of 7 = 10, 20, 50, 100, and 300 
with cold (thin lines) and warm (thick 
lines) electron jets. Figure [Tb shows 
modeled Fermi spectra in vF v units at 
early (a) to late (e) times [Abdo et al. 
2009]. The red lines indicate slope in 
vF v ~ 1 
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sumed to be "Fermi" accelerated at relativistic shocks 
and to have a power-law distribution with index p, 
where N(E e )dE e oc E~' p dE c . Secondly, a fraction £ e 
(generally taken to be < 1) of the electrons are as- 
sumed to be accelerated, and the total electron energy 
is assumed to be a fraction e of the total internal en- 
ergy in the shocked region. Thirdly, the strength of 
the magnetic field in the shocked region is unknown, 
but its energy density (B 2 /8n) is assumed to be a frac- 
tion £b of the internal energy. These assumed "micro- 
physics" parameters, p, £ e , e e and ee, whose values are 
inferred from the spectral fits to the observations [e.g., 
Panaitescu & Kumar 2001; Yost et al. 2003], reflect 
a lack of knowledge of the underlying microphysics 
[Waxman 2006] . It is our intent to place these param- 
eters on a firm physical basis. 



3. SELF-CONSISTENT RADIATION 
CALCULATION FROM PIC SIMULATIONS 

Electrons are accelerated in the electromagnetic 
fields generated by the Weibel and kinetic Kelvin- 
Helmholtz instabilities. Radiation can be calculated 
using the particle trajectories in the self-consistent 
turbulent magnetic fields. This calculation includes 
Jitter radiation [Medvedev 2000, 2006] which is dif- 
ferent from standard synchrotron emission. Radia- 
tion from electrons in our simulations is reported in 
Nishikawa et al. [2011]. 

We have calculated the radiation spectra directly 
from our simulations by integrating the expression for 
the retarded power, derived from Lienard-Wiechert 
potentials for a large number of representative parti- 
cles in the PIC representation of the plasma [Jackson 
1999; Hededal 2005; Nishikawa et al. 2008a,b, 2009b, 
2010, 2011; Martins et al. 2009; Sironi & Spitkovsky 
2009; Frederiksen et al. 2010]. Initially we verified 
the technique by calculating radiation from electrons 
propagating in a uniform parallel magnetic field [Ni- 
hikawa et al. 2009b]. It should be noted that spectra 
obtained from colliding jet simulations (fixed contact 
discontinuity) do not provide spectra in the observer's 
rest frame, and cannot be compared with observed 
spectra [Sironi & Spitkovsky 2009b]. 



The spectra shown in Figure QJ, are for emission 
from jets with Lorentz factors 7 = 10, 20, 50, 100, 
and 300 [Nishikawa et al. 2011, 2012]. In the fig- 
ure we show two spectra for each Lorentz factor (rep- 
resented by the same color line) for initially cold 
(fjet.th = 0.01c) (thin, lower lines) and initially warm 
(wjct,th = 0.1c) (thick, upper lines) jet electrons. Here 
the spectra are calculated for emission along the jet 
axis (9 — 0°). The radiation shows a Bremsstrahlung- 
like spectrum at low frequencies for the eleven cases 
(Hededal 2005) because the magnetic fields generated 
by the Weibel instability are rather weak and jet elec- 
tron acceleration is modest. A low frequency slope of 
vF v — 1 is indicated by the straight red lines. The 
low frequency slopes in our synthetic spectra are very 
similar to those of the spectra in Fig. [Tp in Abdo et al. 
[2009], and show change with the Lorentz factor like 
the temporal evolution observed by Fermi as shown in 
Fig. [Tp. Here we have not included radiation losses 
[Jaroschek et al. 2009; Medvedev & Spitkovsky 2009]. 



4. MAGNETIC FIELD GENERATION AND 
PARTICLE ACCELERATION BY THE 
KKHI 

4.1. Initial Conditions and Previous 
Results 

Recently the kinetic KH instability (KKHI) 
has been investigated using a relativistic counter- 
streaming velocity shear setup with 70 = 3 [Alves et 
al. 2012]. In this RPIC simulation the velocity shear 
occurs at the edges of a velocity field with vq point- 
ing in the positive x (x\) direction in the middle of 
the simulation box, with upper and lower quarters of 
the simulation box containing a velocity field with vq 
pointing in the negative x\ direction as indicated by 
the arrows in Fig. [2f- Initially, the system was charge 
and current neutral. The simulation box dimensions 
were 250 x 80 x 80(c/cj p ) 3 , where w p = {\-Kn\jm, c ) 1 l' 1 
is the plasma frequency, with a resolution of 4 cells 
per c/wp. Periodic boundary conditions were applied 
in all directions. 
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Figure 2: Figure [2^ shows magnetic field lines generated in the relativistic shear from Alves et al. [2012]. Figure [2}} 
shows the electron density in orange (blue) of the plasma that flows in the positive (negative) xi direction. In Fig. [2b 
darker regions in the color map indicate high electron density, whereas lighter regions indicate low electron density. 
Figure 0: shows our simulation model where the sheath plasma can be stationary or moving in the same direction as 
the jet core. 



Figure shows the magnetic field line topology 
for this relativistic case. Alves et al. found that 
KKHI modulations are less noticeable in the rela- 
tivistic regime because they are masked by a strong 
DC magnetic field component (negligible in the sub- 
relativistic regime) with a magnitude greater than 
the AC component,. As the amplitude of the KKHI 
modulations grows the electrons from one flow cross 
the shear-surfaces and enter the counter-streaming 
flow. In their simulations the protons being heavier 
(m p /m c = 1836) are unperturbed. DC current sheets 
which point in the direction of the proton velocity 
form around the shear-surfaces. These DC current 
sheets induce a DC component in the magnetic field 
shown in Fig. The DC magnetic field is dominant 
in the relativistic scenario because a higher DC cur- 
rent is set up by the crossing of electrons with a larger 
initial flow velocity and also because the growth rate 
of the AC dynamics is lower by 7q^ 2 compared with a 
subrelativistic case. It is very important to note that 
this DC magnetic field is not captured in MHD [e.g., 
Zhang et al. 2009) or fluid theories because it results 
from intrinsically kinetic phenomena. Furthermore, 
since the DC field is stronger than the AC field, a 
kinetic treatment is clearly required in order to fully 
capture the field structure generated in unmagnetized 
or weakly magnetized relativistic flows with velocity 
shear. This characteristic field structure will also lead 
to a distinct radiation signature [Sironi & Spitkovsky 
2009b; Martins et al. 2009; Frederiksen et al. 2010; 
Nishikawa et al. 2011, 2012]. 



4.2. Our New Core-Sheath Jet KKHI 
Results 

The simulation setup for our study of velocity 
shears (not counter-streaming shear flows as used 
by Alves et al. [2012] and Liang et al. [2012]) is 
shown in Fig. [2};. In our simulation a relativistic 
jet plasma is surrounded by a sheath plasma. This 
setup is similar to the se tup of our RMHD simula- 
tions [Mizuno et al. I l2007f . In our initial simulation 
the jet core has v COIC — 0.9978c (7 = 15) pointing in 
the positive x direction in the middle of the simulation 
box as in Alves et al. 2012. Unlike Alves et al. the 
upper and lower quarter of the simulation box contain 
a stationary, i> s heath = 0, sheath plasma. Our setup 
allows for motion of the sheath plasma in the positive 
x direction. 

Overall, this structure is similar in spirit, although 
not in scale, to that proposed for active galactic nuclei 
(AGN) relativistic jet cores surrounded by a slower 
moving sheath, and is also relevant to gamma-ray 
burst (GRB) jets. In particular, we note that this 
structure is also rel evant to the jet-in-a-j et or needles 
in a jet scenarios [Giannios et al.l 120091 and papers 
therein], which have been invoked to provide smaller 
scale high speed structures within a much larger more 
slowly moving AGN jet. Similar smaller scale struc- 
tures within GRB jets are also conceivable. 

This more realistic setup is different from the ini- 
tial conditions used by the previou s simu lations with 
counter-steaming flows lAlves et al.l [2012j . and hence 
allows us to compute synthetic spectra in the ob- 
server frame. As mentioned by Alves et al. [2012], in 
a non-counterstreaming or unequal density counter- 
streaming setup the growing kinetic KHI will propa- 
gate with the flow. For GRB jets, the relativistic jet 
core will have much higher density relative to the ex- 
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Figure 3: Magnetic field structure generated by a relativistic electron-ion jet core with 7 = 15 and stationary sheath 
plasma at simulation time t = 70 U)pe- The magnitude of B y is plotted (a) in the y — z plane (jet flow out of the page) 
at the center of the simulation box, x — 500A (d) in the x — z plane (jet flow in the +x-direction indicated by the 
arrow) at the center of the simulation box, y = 100A. Figure [3b shows B y (red), _B X (black), and B z (blue) at x — 500A 
and y = 100A. Figure [3]; shows the x component of the electric current (jet flow is out of the page). The current is 
positive on the core side and negative on the sheath side of the velocity shear. The positive current is stronger than the 
negative current, leading to B y as shown in Figure [3)d. The arrows show the magnetic field in the plane. 



ternal medium. On the other hand, for an AGN jet 
the relativistic core is less dense than the surrounding 
sheath. 

We have performed a simulation using a system 
with (L X ,L V ,L Z ) = (1005A, 205A, 205A) and with an 
ion to electron mass ratio of m\/m c = 20. Figure [3] 
shows the magnetic field structures generated by the 
relativistic electron- ion core with 7 = 15 and with a 
stationary sheath plasma at time t = 70 a;" . Figure 
[3^ shows the magnitude of B y plotted in the y — z 
plane (jet flow is out of the page) at the center of 
the simulation box, x — 500 A. Figure 3d shows the 
magntude of B y in the x — z plane (jet flow in the +x- 
direction indicated by the arrow) at the center of the 
simulation box, y = 100A. Figure 3b shows B y (red), 
B x (black), and B z (blue) magnetic field components 
at x = 500A and y = 100A. Figure [3ji shows the x 
component of the current. Relativistic jet flow is out 
of the page and positive (red) current flows along the 
jet side, whereas negative (blue) current flows along 



the sheath side. Positive currents are stronger than 
the negative currents, leading to B y as shown in Figs. 
[3^ and [3Ji. In forthcoming work we will obtain syn- 
thetic spectra from particles acceler ated by KKHI as 
we have done for shock simulations Nishikawa et al.l 



120111 l2012j| . 



5. Summary and Discussion 

Fig. |4]shows KKHI growth, current filament merger 
and strong magnetic field generation in our simula- 
tion. The structure of KKHI seems similar to that 
found in counter-streaming simulations. In order to 
examine growth rates and wave propagation, we have 
rederived the linear theory to describe our core-sheath 
jets. 

We have extended the analysis presented in Gruzi- 
nov (2008) to core-sheath electron-proton plasma 
flows allowing for different core and sheath densities 
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Figure 4: Figure [4] shows the time evolution of current filaments (J x ) in the area denoted by small box in Fig[3j;. KKH1 
starts to grow at t — 30tJp~ e (a) and current filaments have merged by t = 7Quip e (d). The maximum current density 
(simulation units) is (a) ±0.105 at t = 30Wp e \ (b) ±0.167 at t = 40Wp e \ (c) ±0.296 at t = 50Wp e \ and (d) ±0.625 at 
t = 70Wp~ e . The arrows show the magnetic field {B y , _B Z ) (the length of the arrows are not scaled to the strength of the 
magnetic fields). 



rxjt and n am , respectively, and core and sheath veloc- 
ities zjjt and w am , respectively. In this analysis the 
protons are considered to be free-streaming whereas 
the electron fluid quantities and fields are linearly per- 
turbed. The dispersion relation becomes: 

(k 2 C 2 ± 7a 2 m < am " ^) 1/2 (C - W ara f 

x[(w-fc^t) 2 -^jt] 
± (fc 2 c 2 ± 7 X jt -- 2 ) 1/2 (^-^t) 2 

x[( w -fc^ m ) 2 -< am ]=0, (1) 

where w P jt and CL> pam are the plasma frequencies 
(Wp = 47rne 2 /7 3 m) of jet and ambient electrons, re- 
spectively, k is the wave number parallel to the jet 
flow, and 7jt and 7 am are Lorentz factors of jet and 
ambient electrons, respectively. 



Analytic solutions are not available except in the 
low (ui <C w p and kc <C uj p ) and high (<j 3> w p and 
kc 3> io p ) frequency and wavenumber limits. Equation 
(2) provides an analytic solution to the dispersion re- 
lation in the small wavenumber limit. Here the real 
part gives the phase velocity and the imaginary part 
gives the temporal growth rate and directly shows the 
dependence of the growth rate on the velocity differ- 
ence across the shear surface. 

u ~ (TjtowAyt) , 

(1 ± 7jtWp, am /w p jt) 

± i j^^ /Up f )1/2 M v st - vu). (2) 

(1 ± 7jtWp,am/Wp,jt) 

Equation 2 shows that the wave speed in- 
creases and the temporal growth rate decreases as 
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7jtw p ,am/^pjt = 7jt (^am/^jt) 1/2 increases. Numeri- 
cal solution of the dispersion relation will be used for 
comparison with growth rates and wave propagation 
speeds observed in future simulations. 
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